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Abstract 

Using a self-consistent mean-field theory for the 5 = 1/2 Heisenberg anti- 
ferromagnet Kriiger and Schuck recently derived an analytic expression for the 
dispersion. It is exact in one dimension (d = 1) and agrees well with numerical 
results in d = 2. With an expansion in powers of the inverse coordination num- 
ber 1/Z (Z = 2d) we investigate if this expression can be exact for all d. The 
projection method of Mori-Zwanzig is used for the dynamical spin susceptibility. 
We find that the expression of Kriiger and Schuck deviates in order 1/Z 2 from 
our rigorous result. Our method is generalised to arbitrary spin S and to models 
with easy- axis anisotropy A. It can be systematically improved to higher orders 
in 1/Z. We clarify its relation to the 1/S expansion. 
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The renewed interest in the Heisenberg quantum antiferrromagnet which was raised 
by the discovery of the high T c superconductors and their undoped antiferromagnetic 
mother substances such as La 2 Cu0 4 is still very vivid [1-7] (for a review see ref. [8]). 
There are very few exact facts known about the Heisenberg antiferrromagnet except 
in d = 1 where it can be solved by Bethe ansatz. It is proven, for instance, that the 
translational symmetry is broken in the ground state for all spins S, anisotropics A 
and dimensions d > 2 except for d = 2, S = 1/2 and |A| £ [0.68, 4.55] [9]. 

Besides static quantities such as the staggered magnetization the dynamical quan- 
tities such as the frequency-dependent spin susceptibility x(w,k) are of interest. The 
frequencies at which the latter diverges define the spin wave dispersion relation co(k). 
Furthermore, the transverse susceptibility and the response of neutron scattering ex- 
periments can be deduced from %(w,k). In the traditional analytic approaches for 
calculating %(w,k) which rely on an expansion in 1/5 one does not evaluate k) 
directly but a bosonic Green function. This results either from a Dyson-Maleev trans- 
formation (see e.g. [4]) or from a Holstein-Primakoff transformation (see e.g. [5]). 

Very recently, Kriiger and Schuck proposed a mean-field theory for the spin wave 
dispersion of the S = 1/2 Heisenberg antiferromagnet on hypercubic lattices in d di- 
mensions by using a Dyson equation approach to calculate the dynamical spin suscep- 
tibility [7]. They found a disperson coinciding with the dispersion known from linear 
spin-wave theory modified by a renormalisation of the spin wave velocity Z c given by 
the expression 

Z c\ks = f |f E V 1 " ^( k )) (^ 1-188 for Z = 4). (1) 

V V MBZ / 

Here 7(k) = (1/d) cos(^) and mbz stands for the magnetic Brillouin zone (7(k) > 
0). In d = 1 and for system size iV — > oo (1) yields the correct value of Z c = tt/2. 
For d = 2 the value of (1), Z C \ KS = 1.188, compares very well, too, with other results 
obtained by series expansion (1.18 ± 0.02) [1] or 1/S expansion (1.1765 ± 0.0002) [4] 
and (1.1794) [5] or by quantum Monte-Carlo methods (1.14±0.05) [2] and (1.21 ±0.03) 
[6]. Moreover, we observe that the expansion of the r.h.s. of (1) in powers of the inverse 
coordination number 1/Z, 

Z C \ KS = 1 + ^ + ^-2+°(]p) (^ 1-164 for Z = 4), (2) 

is correct as far as the rigorously known coefficients of the first two terms (orders 1 
and 1/Z) are concerned [8]. 

These facts suggest to pose the question if (1) is rigorously valid in all dimensions. In 
the present Letter we will calculate the coefficient of 1/Z 2 rigorously by an appropriate 
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application of the projection method of Mori-Zwanzig. We will show that the coefficient 
5/8 in (2) is not correct. 

We will use the following appropriately scaled Hamiltonian 

H = jz E (s z a s z b + j(sts b +s a s+)\ (3) 

<a,b> 

on a d-dimensional hypercubic lattice. The sum < a,b > runs over adjacent pairs 
of sites a E A and b G B on the two sublattices A and B. The size of the spins is 
S; the coordination number is Z = 2d; the anisotropy parameter A controls the spin 
fluctuations. We will present our approach for the isotropic case |A| = f to avoid 
that the main ideas are obscured by calculational details, but will state the results for 
general |A| e [0, f]. 

The projection method of Mori and Zwanzig [fO] is designed to compute resol- 
vents, i.e. dynamical quantities, as continued fraction expansions in terms of static 
correlations. We employ essentially the notation of Fulde [11]. An operator set 
{AA] and an antilinear operator product (A\B) are chosen. Here the commutator 
(A\B) := (0|[^4 + , B]\0) with |0) being the ground state is the appropriate choice [12]. 
Note that this is not a positive semi-definite scalar product. The minimum operator set 
we use is A = {A 1 , A 2 } with the spin-flip operators A^k) = sJl/N Y,i^A S~ exp(ikR,) 
and A 2 (k) = ft/N^g S~ exp(ikR i ). 

The Mori-Zwanzig method is applied to the time-dependent spin-spin correlation 
function Rij(t) := — i(T{ Af(t)Aj}), where T is the time ordering operator. The 
causal 2x2 matrix Green function Rij(t) contains all information about %(w,k). 
With the Liouville operator L := [H, ] the Fourier transform Rij(co) equals Rij(co) = 
(Ai [uj — Z/p 1 Aj^j. The resolvent R(w) can be rewritten [11] 

R(w) = — (L + M(w))P _1 ] 1 P , (4) 

where bold face letters stand for matrices. The matrix elements of P,L and M(w) are 
defined by 

Pij := (Ai\Aj), L id := (Ai\LAj), M^(uj) := (Q A LAi \[co - QaLQa]' 1 QaLAj) , (5) 

where Qa is the projection operator onto the operator subspace orthogonal to A = 
{A 1 ,A 2 }, i.e. Qa-=1-T, \Ai) (P" 1 )^ (AA. 

One obtains an approximation for R(w) by neglecting the memory matrix M(w). 
This approximation is improved by including more and more operators in the operator 
set considered [11]. The basic Lanczos idea that new operators are generated by the 
iterated application of the Liouville operator L serves as a first guideline to which 
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operators should be included. The spin wave dispersion co(k) is determined from the 
poles of the approximate expression for R(w). 

The memory matrix M(w) is of the same form as R(w). With A\ := Q A LA { and 
V := Q A LQ A one can define matrices P', L' and M'(w) analogously to (5) and express 
M(ci>) according to (4). If we are able to calculate the order in 1/Z of the memory 
matrix M(w) we know in which order an approximate expression for R(w) neglecting 
M(k>) is exact. In practice, we will compute the iterated projector matrix P' of M with 
P[ - := (QaLA^QaLAj) and determine its order in 1/Z. This order obviously is the 
order of M(w) for large oo. It is also the order of the error of the spin wave dispersion 
u)(k), which is determined from the approximation for R(w), because M(w) does not 
have poles around co = oo(k) [4]. 

To illustrate the procedure described above we calculate for the set A 

where u := — (2{S^S b ) + (S+S b ~))/(2S) and a £ A and b G B represent adjacent 
sites on the two sublattices. The 1/Z expansion of the expectation values results in 
(S*) = S - 1/(4Z) - 9/(16Z 2 ), v = S + (3 - 6S)/(16SZ 2 ) and ac = -(1 - 7 2 (k))/Z 2 . 
From ac = (9(1/Z 2 ) we know that we may trust the dispersion relation resulting from 
the approximation for R(w) based on {^i,^} only to order 1/Z. In this order one 
obtains co(k) = (1 + l/(4S , Z)) v /l - 7 2 (k). 

To derive (1), Kriiger and Schuck inserted (6) into (4), neglected M(cti), and used 
certain self-consistencies specific for S = 1/2 to determine the expectation values v 
and {S*) [7]. Thus one understands why (1) yields the coefficients up to order 1/Z 
correctly. 

The 1 /Z expansions of the static expectation values needed in the projection method 
can be computed by ordinary perturbation theory in the spin fluctuations starting from 
one of the Neel states. To be exact up to order 1/Z m one calculates the perturbation 
series in A to order A 2 ™. Furthermore, the 1/Z expansion simplifies complicated 
integrals over the Brillouin-zone. Any expression containing 7 is expanded in powers 
of 7 which can then be easily summed. In the present calculations we performed these 
expansions using an algebraic computer programme (MAPLE). 

In an attempt to describe the dynamics properly up to order 1 /Z 2 what comes into 
mind first is the simple Lanczos extension A' = {A 1 , A 2 , A[, A' 2 } of the operator set 
{A,,A 2 } with A\ = CE W) e' kR «(5 6 5: - S;S§) and A' 2 = CJ2^ b) e^(-S b S z a + 
S~S b ), where C := Qa\J2/N j{SZ). Inspection reveals, however, that the correspond- 
ing memory matrix is still of order 1 / Z 2 and not 0(1/ Z 3 ) as desired. Thus this straight- 
forward extension is not sufficient. 



L A = 2u 



1 - 7 (k) 

- 7 (k) 1 



P' A = K 



1 
-1 



(6) 
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Let us consider a general operator set K. To describe the dynamics in order 1/Z 2 
correctly it is necessary that for any operator O E K the deviation of LO from K is 
0(1/Z 3 ). Splitting the Hamiltonian H = H MF + Hi in H MF := - Y, a £A S z a + E&es S z b + 
NS/2 and H, := (SZ)- 1 E< a , fc > [(S z a - S)(S* b + S) + (A/2)(S+S b ~ + allows to 

distinghuish between the part of order 1 (H MF ) and the part of order 1/Z (H^ since 
the ground state approaches the Neel state on Z — > oo. Commutation with H MF does 
not increase the order in powers of 1/Z. Therefore, we require that K is closed under 
the Liouville operator L MF = [H MF , ]. Inspecting A\ and A' 2 one realizes that they are 
not reproduced by L M f because the relative signs of the two terms on the r.h.s. are 
changed. This leads us to the set B = {B 1 , . . . , B 6 } with B 1 := A 1 , B 2 := A 2 and 



B 3 :-- 
B 5 :-- 



«kR„ 



<»,*) 

-c e s b s:^ 

(a,b) 



B A 
Bq 



C E S a S z b e^ 

(a,b) 

C Y: S b S z a e> kR ° 

(a,b) 



(7a) 
(7b) 



The operators B 3 , . . . , B 6 are the split parts of A[ and A' 2 : A\ = B 3 + B 6 and A' 2 
B4 + B$. They are eigenoperators of L M f, L MF Bi = ±Bi, whence Qb^mp^ = 0. 
Explicit calculation to order 1/Z 2 yields for the operator set B 
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N : = 



1 + 7 2 

-27 



-27 
2 



, (8) 



where 7 is short for 7(k) and the matrix V := —[[0, 1], [1, 0]]. We find a block structure 
with 2x2 blocks. The upper left block represents the information that we had already 
with the first set {A 1 , A 2 }. The matrix elements of P# which link B 1 , B 2 with the other 
operators are zero by construction since QaB\/2 = 0. The corresponding lateral matrix 
elements of L B are easily obtained from (B^LB^ = (Bi\B 3 + B 6 ) and (i^|L52) = 
(Bi\B4 + 5 5 ) for i > 2. We expect the simple block structure of the lower right 4x4- 
block to disappear in order 1/Z 3 . 

Explicit evaluation aided by computer algebra shows that V B of M B is 0(1/ Z 3 ). 
By the line of argument presented above we know therefore that the dynamics is now 
correctly described to order 1/Z 2 . 

Using (8) and neglecting M B (u) we obtain a 6 x 6-resolvent. From this we use 
the upper left 2 x 2-block as improved approximation for R(w) which is correct to 
order 1/Z 2 . The zeros of det — L^P^ 1 ] give the dispersion relation co(k). One finds 

Z 2 C (1 - 7 2 (k)) + 0{1/Z 3 )) with 



det [to - L^P^ 1 ! = 



u 2 - 9 + 0(1 /Z)f(u 2 



Z c = 1 + 



+ 



ASZ 16SZ 2 



+ O (-^j (« 1.148 for Z = 4, S = 1/2). 



(9) 
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With this renormalisation factor co(k) = Z c yjl — 7 2 (k) represents the systematic sec- 
ond order 1/Z expansion of the spin wave dispersion. To this order the spin fluctuations 
only renormalise the spin wave velocity; they do not affect the shape of the dispersion. 
We do, however, expect a modification of the shape from higher orders of the 1/Z 
expansion. 

For S = 1/2, (9) differs from the Kriiger/Schuck result in (2). We therefore come 
to the conclusion that the analytic expression (1) cannot be rigorously valid in all 
dimensions, not excluding completely the possible validity in any particular dimension. 
There is no fundamental reason known why the mean-field theory of Kriiger/Schuck 
yields the exact results in d = 1 although a similar agreement was already found before 
[13]. 

Investigating the pole structure given by det — L^P^ 1 ] one finds that besides 
the dispersion pole there is a double pole at u = 3 which is dispersionless. Physically 
it corresponds to 3-magnon processes. They are dispersionless since we desribe the 
3-magnon-pole only in leading order where due to ^1 — 7 2 = 1 + G(1/Z) on averaging 
in the Brillouin zone the magnon energies are replaced by 1. 

It is obvious that the approach described above can be extended to the model (3) 
with an easy- axis anisotropy |A| < 1. We state here the final second order result for 
the dispersion relation 

^ (k) = ,_ AV+ A ! (' + AV-2f) + 3A'(3AV-4 T ' + l) 



2SZ ' 8SZ 2 

A 2 (9A 6 7 4 - 3A 4 7 2 (4 7 2 - 29) + 4A 2 ( 7 4 - 52 7 2 + 6) + 128 7 2 - 32) / 1 



+ y .do) 



16S 2 Z 2 (8 + A 2 7 2 ) 

There is an interesting relation between the 1/S expansion and the 1/Z expansion of 
the spin wave dispersion: the correction of order S~ m in the 1/S expansion calculated 
for a hypercubic lattice of arbitrary dimension is small of order Z~ m . This implies that 
taking this 1/S expansion up to the order 5""™ one obtains all contributions up to the 
order Z~ m of the 1/Z expansion. This also explains why in our 1/Z expansion terms 
of order Z~ m contain only contributions of order S~ k with k < m. 

Note that the hydrodynamic relation Z 2 = Z p /Z x [8] is intrinsically fulfilled since 
all three renormalisations can be derived from R(w). With the dynamic transverse 
susceptibility k) := — + 2i?i i2 + i?2,2)/4 one finds 

Z x = limk^oli m w^o(l/Xciass)x(w, k) and for the spin stiffness renormalisation Z p = 
- lim^o lim k ^ (l/p class )(w 2 /k 2 )x(w, k) [14]. Our results are Z x = 1 -1/(2SZ) + (5 - 
185) /(2AS 2 Z 2 ) + O (Z- 3 ) and Z p = 1 + (1 - 185) /(48S 2 Z 2 ) + 0(Z~ 3 ). 

In conclusion, a systematic 1/Z expansion of the dynamical spin correlation func- 
tion R(w) up to the order 1/Z 2 was presented. Explicit results for arbitrary spin S 
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and easy-axis anisotropies |A| < 1 for the spin wave dispersion up to the order 1/Z 2 
were given. The result for S = 1/2 and |A| = 1 was used to conclude that the spin 
wave velocity (1) of Kriiger and Schuck is not exact in all dimensions. A general recipe 
how to extend the 1/Z expansion to higher orders was given: One has to make sure 
that the set of operators generated by iterated application of the Liouville operator L 
is closed under L MF . Finally, the relation of the 1/Z expansion to the 1/S expansion 
was elucidated. 
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